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ABSTRACT 
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results are compared with other related publications. 



*Present address: High Energy Physics Institute, Tbilisi State University, University St. 9, 380086 
Tbilisi, Republic of Georgia. 

^e-mail: mereb@sun20.hepi.edu.ge 

^e-mail: apcont@physics.mcgill.ca, acontog@cc.uoa.gr 

•^e-mail: ggrispos@cc.uoa.gr 



I. INTRODUCTION 



Deep inelastic scattering of longitudinally polarized particles has provided important 
information on the spin structure of the nucleon. However, the size and shape of the polar- 
ized gluon distribution At? in the proton remains an essential problem. Significant progress 
requires experiments on reactions with longitudinally polarized particles dominated by sub- 
processes with initial gluons. Such a reaction is 

7 + P^Q(Q) + x, (1.1) 

where Q{Q) denotes heavy quark (antiquark); this is dominated by 

1 + 9^Q + Q (1.2) 

An experiment closely related to ( |l.lj ) is soon going to take place [0] and there are more 
than one proposals 0. 

At the Born level, ( p. . 1| ) has been studied long ago 0,0. However, the importance of 
knowing the next-to-leading order corrections (NLOC) cannot be overemphasized. This 
work presents detailed results on a NLOC calculation. 

It should be noted that NLOC for ( |1.1| ) have already been published 0. We believe, 
however, that in view of the importance of (1.1), an independent determination of NLOC 
in a different regularization approach (see below) is in order. Extensive comparisons with 
the calculation of |^, as well as certain differences in our view regarding certain questions 
will be also reported. 

At NLO, apart from the loop and gluon Bremsstrahlung (Brems) contributions to (|1.2| ), 
the subprocesses 

i + q{q)--Q + Q + q{q), (i-3) 

where q denotes a light quark, should also be taken into account. 

We note that the Abelian part of NLOC for (1.4) provides the corrections to 

1 + + Q (1.4) 

This part has already been determined NLOC to (1.4) are of interest in themselves 

in connection with Higgs boson searches when the Higgs mass is in the range of 90 to 160 
GeV. 

The loop and 2-^3 parton graphs involved in NLOC introduce ultraviolet (UV), 
infrared (IR) and collinear singularities, which are eliminated by working in n = 4 — 2e 
dimensions. For polarized reactions this requires extension of the Dirac matrix 75 in n 7^ 4 
dimensions. Unless otherwise stated, we work in the scheme of dimensional reduction 
(RD), which simplifies the calculation of the traces. Certain subtleties of RD have been 
discussed in and are mentioned below. Furthermore, we use parton distributions whose 
evolution, via 2- loop anomalous dimensions, is determined in a scheme different from RD. 
This necessitates the addition to our perturbative results of certain conversion terms. 
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In all the above contributions the photon interacts in a direct way. In addition, there 
are also resolved contributions, in which it interacts through its partonic constituents; in 
fact, strictly speaking, at NLO, scheme independent cross sections arise only by adding 
them. At this moment a complete calculation of the resolved contributions is not possible, 
and we will be limited in giving an estimate. 

The paper is organized as follows. Sect. II contains our general procedures. Sect. Ill dis- 
cusses the loop contributions to the photon-gluon fusion subprocess and Sect. IV the corre- 
sponding Brems ones. Sect. V presents analytic results on the subprocess (1.3). In Sect. VI 
we derive the nesessary formulas for calculating various physical observables. Sect. VII 
presents our numerical results and discusses the possibility to distinguish between three 
sets differing essentially in the polarized gluon distribution function Ag. Sect. VIII deals 
with our comparison with [^], as well as with Sect. IX presents our conclusions. Finally, 
in three Appendices we present results completing our determination of NLOC 



II. GENERAL PROCEDURES 

The Born and the loop contributions to 'j + g^Q + Q are shown in Fig. 1. With the 
4-momenta pi,i = 1, ...,4, as indicated and with m the heavy quark mass we define: 

s = {pi+p2)'^, t = T -m^ = {pi-p^Y -m^, u = U - = {p2 - P'if - (2.1) 

Let Mj(Ai, A2) the amplitude of any of the contributing graphs, where Ai, A2 the helicities 
of the initial partons; our polarized cross sections correspond to the quantities: 

\nM^{++) 4- (++) - M, (+-)] (2.2) 

where S denotes summation over the helicities and colors of the final particles and average 
over the colors of the initial. For the determination of the asymmetries we need also the 

unpolarized cross sections, which correspond to the average of Mj(++) Mj (++) and 

M,(+-) M, (+-). 
We also introduce 

v = l + t/s w = -u/{s + t) (2.3) 

To reduce the length of the subsequent expressions we will make use of the results presented 
in 1^. Thus our leading-order (LO) polarized and unpolarized cross sections are 

^ ^ dvdw ^ dvdw ^ ^ 

where k = as/Saeg and [A]dai^o/dvdw the corresponding [polarized] unpolarized cross 
sections for 77 — > QQ (Eq. (9) of p). For later use we note that [A]d(7i^o/dvdw are 
proportional to: 



AB{s,t,u) 



tu 
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and (see also ||^) 

B{s, t,u) = -( — + 4^ 1 - ^ (2.5) 

s \ tu tu \ J J 

In determining the loop contributions, the renormalization of the heavy quark mass and 
wave function were carried on shell, as in P], i.e. the renormalized heavy quark self-energy 
Er(p) was taken to satisfy at = rn?\ 

K{p) = ^^-^P^ = (2.6) 

This determines the mass and wave function renormalization constants Zm and Z2 [0. 
Dimensional reduction does not automatically satisfy the Ward identity 

Zi = Z2, 

where Zi is the renormalization constant for the vertex of the graph Fig. 1(d). This requires 
the introduction of proper finite counterterm, of which the form is given in 0. 
In the present case charge renormalization is also required. Defining 

let go{g) be the bare (renormalized) coupling, Zg = qq/ g the charge renormalization con- 
stant and b = (llA'^c — 2Nif)/Q, where Nif is the number of light flavors. We take 

Zg = l-j [c,iM)b - ^C,(m)| (2.8) 

where M is a regularization mass. In this scheme the contribution of a heavy quark loop 
in the gluon self-energy is subtracted out, i.e. the heavy quark is decoupled This 
is consistent with parton distributions AFa/p{x,Q'^) of which the evolution is determined 
from split functions involving only light quarks, as is the case of AFa/p used below. 

Finally, the renormalization of the gQQ vertex was carried using the Slavnov- Taylor 
identities iTDI . 



III. LOOP CONTRIBUTIONS 



The loop graphs contributing to (1.2) are depicted in Fig. 1. The integrals for the 
Abehan type of graphs (a)-(e) were calculated in 0. The non-Abehan graphs (g) and (h) 
introduce tensor integrals of the form 



(27r)" q^{q - P2Y[{q + Pa- P2Y - m 
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and 



q^^, q^q'^, q^q'^qP 



As in Q, using Passarino-Veltman techniques |]rT|, we reduce them to scalar ones; those 
can be found in [0. 

The contributions presented below include the t ^ u crossing symmetric of Fig. 1 plus 
UV counterterms; thus they contain no UV singularities. 
The graphs (a)-(e) give: 



Nc daa- 



dvdw dvdw 2 

where damsel dvdw is given in Eq.(16) of and 



dv 



'■-5{l-w) 



(3.1) 



dan 



dv 



ira2Ai[2(C(2)-Li2(— ))(l + 3 



t 



ln( 



;)il + —) + 2]+A,ln{- 



T 



+yl3(Li2(^) - C(2)) + + (t ^ u)] 



(3.2) 



with 



Kr = —aa^e 



2„2 



8s ^ 

Here and subsequently the polarized cross sections are given by (3.1) and (3.2) with 
Ada/dvdw, Ada/dv and AAi, i=l,...,4, replacing the corresponding unpolarized quan- 
tities. The [A]Ai are given in Appendix A. 
Graph (f) contributes: 

daj ,^ Nc^dabox 



dvdw ^ ^ 
with da^ox/dvdw in Eq.(22) of 0. 

Turning to the non-Abelian graphs, (g) gives 



2 dvdw 



(3.3) 



where 



with 



and 



dal_ 
dvdw 



dag 
dv 



-2TxasCe{m)N^ 



c- 



dvdw \e'^ e 



Nc ( ddn 



+ 



dan 



2F{e)\~A,-\ni^) + A\ 



dv dv 
t 



6{l-w), 



e"^ 6 6 IV?' 



u 



tu — sm? 



da, 



T 



-t 



T 



Y = K,{2A,[U,{ — ) + ln^(— ) - 2] + A\[AU,{—) + 41n^( — )] 

+A'2H^) + A'^ + {t^u)} 



(3.4) 

(3.5) 
(3.6) 



(3.7) 



In (3.5) and (3.7), [A]Ai are given in App. B of and [A]A'i, i=l,2,3, in App. A of this 
paper. 

The contribution of the graph (h) is: 



dvdw 



c 



\-nci C (m) — ^^^ ^ fdoh^doh 
^ ^ dvdw \ dv dv 



5{l-w) 



(31 



where 

doh _ 
dv 

and 



2F(e) 



ln(4) + 21n(- 



1 2 2, -t. 



+ {t^u) \ (3.9) 



^ = i^UA.hf C(2) - Li,(^) + 41n(4) ln(4) - ln^(4; 



-t. 



-t. 



+ B\U2{ 
-t. 



+ {B'2 + ^'i)C(2) + B', ln^( — ) + BU ln(— ) + B', ln(— ) ln(— ) 



-B'e + {t^u)} (3.10) 



The coefficients [A]B'i, i=l,...,6, are given in App. A. 

Finally, after cancellation of the UV singularities, graph (i) does not contribute. 

We remark that regarding the terms the contributions of the graphs (g) and (h) 
taken separately are not proportional to the Born dajj^/ dvdw\ only their sum is propor- 
tional to the Born. The same holds regarding the terms 1/e. 



IV. GLUON BREMS CONTRIBUTIONS 



In this chapter we present complete analytic results for the NLOC arising from Brems. 
To the best of our knowledge, in relation with heavy quark production, such results have 
not so far been presented. 

With k the 4-momentum of the emitted gluon we introduce also 

S2 = {k+ p^Y — = s + t + u = sv(l — w) (4.1) 

The Brems graphs contributing to the NLOC of (1.2) are shown in Fig. 2A. The squared 
sum of the corresponding amplitudes (plus those obtained via pi P2) after summing over 
final spins and colors and averaging over initial colors is given by 

Am^mU' = Ks{e) {^G^^ - ^G^^) (4.2) 
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where G^'^ is the quantity in the square bracket of Eq.(24) of (plus pi <-> P2), G^^ has 
the expansion 

Cx^y = d H ^ 5- H ^ H H h 68 

P2 - Pa Pi- Pi Pi ■ Pi Ps ■ k Pi ■ Pa Pz ■ k pi ■ Pa P2 ■ Pi Pi- Pi 

eg , P2-k ^ en , . Ps ■ k . ps ■ k"^ ~ ■ k'^ ~ ■ k 



+ r + eio^ + r — r + fi- — + f2- + /3^^^ + /. 



P2-PiP3-k ' '^'"ps-k ' pi-pip3-k ' ''^P2-Pi ' ''^P2-Pi ' ''^ P2-k ' ''^p2-k 
_^ /s _| /e |_ fi _| is |_ /g , /lo 



j02 ■ A; Pl ■ Pi P2- k P2- k"^ Pl ■ Pi P2 ■ k"^ P2 ■ k ps ■ k Pi ■ pi P2 ■ k 

+ i^^^ (4-3) 

Pl ■ P4 P2 - k^ 

and 

As in Sect. Ill, AjMs^lgp is given by (4.2) and (4.3) with AG^^ , AG^s, Ae^ and Afi, 
i=l,2,...,13, replacing the corresponding unpolarized quantities. The coefficients [Aje,, 
[A]/j of (4.3) are given in App. B. 

The Brems contribution to [A]dcr/dvdw is obtained by working in the Gottfried- Jackson 
frame of Q{Q) and gluon (cm. system of p^ and k). Details are given in 0. The terms 
with coefficients [A]ej, [A]/j in (4.3) give contributions singular at S2 = (w = 1) and 
must be integrated in ri 7^ 4 dimensions. In view of the fact that the 2^3 particle phase 
space IS proportional to sl~^^ (Eq. (26) of 0), the remaining terms can be integrated in 
4 dimensions. The arising integrals are given in [|l^]. Certain terms of special interest not 
given in |^ are determined in App. C. 

Corresponding to the second term in (4.2), with y = \J (t + — Am?s^ 82 = 82 + m? 

and X = (1 — /3)/(l + /5), where (3 = y^l — Am'^/s, the final result is 
d< Kb{0)Nc^ VS2 . 282,82 , 4^2 , 252 , ^2 

+ 67/8 + eg/io + eio/i6 + eii/i3(t ^ u) + fiFi + /2F2} 

-F(0)|— In 65 + eehiit ^ u) + eg/nj 



(1 -w7)+52 16 '^52^ T + U + y 

- ^t|^(0){/3^3 + + hF,' + fsF^ + fioF^o + /ii^n - (2 In ^ + In ^)(/3F3^ 

O2 iO TTl 02 

+ fiFl + fsF^ + hoF^, + fnF^,)} 

- ri if + hF^ - (2 In ^ + In '^){hF^ + hF^ + hF^)} 

(1 — w)j^ D2 16 D2 

+ 2L^'l^FmhF^ + UF^ + ~hFI} 
k>2 16 
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+ 87ra,iVcC,(m)[A]^{21n2 ^ - ln\x) + lln\-x) + 21n^ln^ + U,{1 - 

avaw m"^ 2 u t x t 

In Eq. (4.5), the integrals /j are given in the App. C of and the integrals Fi in the App. 
C of this paper. Also, L+ = (ln(l — w)/{l — w))+, which enters through the relation 

(1 - wr'~'^ = -^5(1 -w) + J^^^ - 2^^+, (4.6) 
where the so called "plus" distributions are defined in a usual way: 

■ ,,fM_ ^ r (4.7) 

(1 — w)+ Jo 1 — w 

For the second term of (4.2), we give the terms ~ and 1/e, as well: 



-rt = ^8™sC',(m) A -^-^87ra,C',(m) A -^{31n(— )-ln(— 
avaw avaw e avaw 



+ ^87ra,C,(m) A ^ln(a;) 

1 2s7; 

- -^3^F(e)[A]P/^(x2)[A]S(x2.,t,X2«) (4.8) 

where [A]Pg^g(a;) is the 4-dimensional g gg split function without the 5{l — w) part, -F(£), 
[A]i?(s, t,M) are given by (3.6) and (2.5) and 

X2 = (4.9) 
1 — vw 

Addition of loop and Brems contributions cancels the singularities and part of the 
1/e. The remaining 1/e are cancelled by a factorization counterterm corresponding to the 
final gluon emitted coUinearly with the initial one (Fig. 2 A, graph (d)). In the MS scheme 
this counterterm gives: 

= i^f(e)|A]P„(..,|A|B(..M...«) (^J, (4.10) 

[A]Pgg{x) the g gg split function and Mp the factorization scale. 

Our cross sections will be convoluted with parton distributions evolved via two-loop 
split functions. In n dimensions the split functions have the form 

[A]A"„(a:,£) = \/\]PUx) + s[/\]PUx) (4.11) 
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The polarized split functions have been determined |T3| , p!^ in the t'Hooft-Veltman scheme 
[|15[] modified so that AP^(a;, e) = Pqq{x,e). In this scheme 

AP;^(x) = 4iVc(l -x) + - x) (4.12) 

However, our calculations were carried in dimensional reduction (Sect. I), where 

^Pl,{x) = P:,{x) = (4.13) 

Thus, a conversion term Adaconv/dvdw should be added to our Ada^^ /dvdw. Conversion 
terms are determined from the difference of [A]P^^(x) in the two schemes |jl^: In the 
present case 



^ conv ^ _ F(0) APyXs) AP(X2S, t, X2U) (4.14) 

dvdw 1 — vw 
with APgg{x) given by (4.11). 

The unpolarized parton distributions we use were evolved in the MS scheme, where 
Pgg{x) = 0. Thus conversion term is not required. 



V. SUBPROCESS QQq 

The graphs contributing to this subprocess are shown in Fig. 2B. The squared sum 
of the corresponding amplitudes, after summing over spins and colors and averaging over 
initial colors, is given by 

Am'\M]l,\' = -^(ATifaalielQ, + ejQs + ege.Qs) (5.1) 

where the charge of the light quark q. |M2^3p corresponding to 7g QQq is given by 
the same expression with an opposite sign of the last term. The quantity Qi is given by 
an expansion similar to (4.3). Next we introduce 

-534 =P3-p4 + m^ (5.2) 

Then Q2 and Q3 are of the form: 

Q2,3 = ei H h eg h /4 r H r H r + fi2 h /i3^— 

Pi - Pi Pi - Pa P2-k P2-k Pi-PiP2-k S34 s|4 

, /l4 /l5 , /16 , /l7 , /18 , /l9 . P P3 ■ k 

+ ' T T + ^ T ^ 1 T + .^20 r 

S34 Pi ■ S34 pi ■ fc s^4 pi ■ fc S34 Pi ■ P4 pi-k p2-k Pi - k 

+ — (5.3) 

S34 P2 ■ 
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As before, A|M2f,3p is given by (5.1) and (5.3) with AQ^, r=l,2,3, Acj and Afj replacing 
Qr, Cj and fj. The coefficients [Ajcj and [A]fj are given in the last part of App. B. 

The contribution to [A]da/dvdw is obtained by working as in Sect. 4 (cm. frame of 
Q{Q) and final light quark). Again the terms with coefficients [A]ej and [A\fj must be 
integrated in n dimensions. 

After phase space integrations, we get the following results for the sets Qi, Q2 and Q3: 

= Le|{ei + ^-^— 63 + —^ In + 68/10 + /4i^^ + Zeiss' 

dvdw ^ m"=(s + rj^ S2[s + t) m"^ 

+ f,F^ + f,,F^^ + f,,F^^} (5.4) 



,Q2 



with 



= Le2{ei + /12F12 + /13F13 + /14F14 + fwF^e + fnFii + /2oi^2o} (5-5) 

^ - ^e,e,{e. + ^ln-Ae4 + e3/.o + /.Fr + /eF,^ + A2F. 

+ /i4i^i4 + /lei^^e + /isi^is + /igi^ig + /2oi^20 + /2i^2i} (5-6) 

2 1 S2 



' Nc Stt ^2 

We do not write down expressions containing 1/e poles coming from sets Qi and Q2 as 
they are equal with an opposite sign to the corresponding counterterms with (S'2m^/s|)^ 
instead of {w? /MpY (see below). 

The singularities arise when the final light quark is collinear with the initial one {k-p2 = 
0, Fig. 2B, graphs (a), (b)) as well as when the photon is collinear with the light quark 
{k ■ pi = 0, Fig. 2B, graphs (c), (d)). To eliminate them we introduce two counterterms. 
In the second case the counterterm involves the Born cross section for QQ, which is 

proportional to p!7|: 



ABgg{s, t, u) = -B,g{s, t, u) = ' ' [ + 2" l (5.7) 

Moreover, in n dimensions, in the t'Hooft-Veltman scheme: 

AP;^(x,£) = Cp{2 -x + 26(1 - x)} P;,(x,£) = {i±ii:i^ - ex} 

and 

AP;^(x, e) = x-^-Eil-x) Pj;(x, e) = + i^^- _ ^^(^ _ ^) 



Thus, in the MS scheme, the ffist counterterm gives: 

r . T dail^ 1 2SV , X r . n / N r . n / N / \ ^ 



where F{e), X2 and [A]B given by (3.6), (4.8) and (2.5), and the second gives: 

'^'^ = ~m§JA]P„MlA]B,,i^s,^t.u) [^J (5.9) 

Although not necessary, it is now customary and even advantageous ^ to average the 
unpolarized cross section over n — 2 spin degrees of freedom for every incoming boson. This 
convention is employed when fitting the unpolarized structure functions. As a result, for 
the unpolarized case, the r.h. side of (5.8) should be multiplied by (1 + e) and of (5.9) by 

(1-5). 

We note that, upon integration, the singular terms in [A]Q3 cancel out, as they should 
since there is no counterterm proportional to CgCQ. 

Conversion terms are also needed in the present case. Along the lines of the previous 
section: 

r/rr^^) 2 s?; 

[A]-^ = --^F{0)[A]P^^{x,)[A]B{x,s,t,x,u) (5.10) 
avaw 1 ~ vw 

and 

= -^mf[m,{w)[A]B,,{ws,wt,u) (5.11) 



Finally, we have carried our analytical calculations using REDUCE |T^ and to some 



extent FORM 19 



VI. PHYSICAL CROSS SECTIONS 



Here we present the necessary formulas needed for calculation of the differential and 
total cross sections for the physical process jp Q + X . This includes derivation of 
physical cross sections for both components of the reaction, i.e. pointlike and resolved, 
the latter to leading order (LO). Note we always observe a heavy quark in the final state. 
Capital letters in this chapter refer to the kinematic variables of the physical process and 
small letters to those of the subprocess. Starting with the pointlike component, the total 
cross section for the reaction (IJ.) can be written 



[A]a,p{S)=f' dx[A]U/p{x,Q')[A]a^,{sy, (6.1) 

b denotes the corresponding parton, f^/p its probability distribution and 

s = xS, Xmin = 4:rn'^/S. (6.2) 
The total partonic cross section [A](3"^fc(s) can be calculated straightforwardly by 

|Al*,.(,)=r"A,/' My,^^-) ^ (6,3) 
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where 

1 771^ 

'^max/min T^i.^ Wmin JZ T- (6-4) 

2 sv(\ — V) 

To derive the transverse momentum differential cross section we note that the transverse 
momentum pt of a heavy quark is invariant under boosts along the beam axis; also, that 
our heavy quark rapidities are defined with respect to the photon. Consequently we have: 

da^^^r^ |A1A,,(..Q=)|A|;?^2^. (6.5) 

with 

XminiPr) = KPt + (6.6) 
da^b{s,pT) ^ r?^— 2pT da^b{s,v,w) 
dpT Jymin sv dvdw 
The integration limits on the cm. rapidity y are 



Vmax = -Vmin = + J wj - 1 



with Wm = 4(p|. + m )/s. Integration over rapidities in (|6.7|) is not well defined for "plus" 



distributions (given in ([4.7| )) in the partonic cross section. The problem is solved with a 
change of variables; One needs to consider an integration contour for a heavy quark rapidity 
y and split it into two parts that no overlapping (i.e. double counting) occurs. Formally 
one would have 

dyfiy) = / dyfiy ^ -y) + dyf{y). (6.8) 

In particular, the splitting point yo must be the point where function w = f{y) has a 
minimum. We find 

yo = \n{2^J{pl + m^)/s) 
and the general relation between the "new" variable w and the old variable y is 



The correct sign in ( |6.9[ ) is different in different integration regions, e.g. in the region 
[yo,ymax] the function e~^ decreases when one goes from to ymax, thus minus sign in 
( |6.9|) . Similarly, we find that for the first term of (|6.8|) the sign for e~^ in (|6.9|) should be 
positive. The resulting expression for the pt differential cross section reads: 

da^b{s,pT) dw 2pT f, . da da \ 

A — = / ; <^ [A]——{v = v+) + [A]——[v = V-)\ , 6.10 

apT Jwm w g h _ Wm/w [ dvdw dvdw J 
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where 

v^ = ]^{l±^l-wjx). (6.11) 

However, even the expression ( |6.10|) is not well suited for numerical integration. One 
notices that there is a numerically divergent (though analytically integrable) square root in 
the denominator. The singularity comes from the lower limit XminiPx) of the x integration. 
To avoid this minor problem one more change of variables is necessary. Instead of the old 
variables x, w we introduce the new variables z, w' through the relations 



V wx 



w' = w. (6.12) 



To correctly define integration limits for the new variables one has to perform a nontrivial 
mapping. Finally we obtain: 

(6.13) 

Fb/pix^Q"^) is a momentum distribution and 

Zm = ^l- XminiPT), = — X = . (6.14) 

For the rapidity Y fixed one gets the following expression: 

da^^jS^ = /' dx [A]/,/,(x,Q^)[A]^^^^^, (6.15) 
dY JxmUY) dy 

x„.^n{Y) = e-^/(v^/m - e^), (6.16) 

da^b{,s,y) ^ /•! 2wdw da 
dy Jwrnin (e^ + we~yY dvdw 

^sjm - e-f ' 1 + WQ-^y ' ^ 2 

Finally we turn to the resolved LO photon contributions. We define the doubly differ- 
ential cross section da/dYdpx for the 2^2 subprocess: 

rfo^ ^ , [A]Fa/^(xi, Q^) [A]Fb/p{xl Q^) r ^ 1 . . . on 

A f = 2pT / dxi '-^ [A](Ta6(s,Xi,x^), (6.18) 

where 

x^^^i^, A^{^±^Y/\ s = x,xlS. (6.19) 
Xi — e /i o 

The expressions (6.1) - (6.19) give all the formulas we have used. 



with 
and 
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VII. NUMERICAL RESULTS 



We present results for Q=c-quark {rric = 1.5 GeV) at ^S^p = ^/S = 10 GeV, relevant 

to the experiments |1]| and [0, (a)] and \^ = 100 GeV, as well as for Q=b-quark (m;, = 5 
GeV) at 75 = 100 GeV; the later energy is relevant to HERA. Higher HERA energies are 
not considered as the cross sections become too small. The effect of changing rric is also 
considered. 



We use the NLO sets of polarized parton distributions of |2y], which can be characterized 



in terms of the polarized gluon distribution Ag ( follows: 
Set A: A^f (x) > and relatively large 
Set B: A^f (x) > and small 

Set C: Ag (x) changing sign; Ag (x) < for x > 0.1. 

Notice that in the presented results also the LO contribution is convoluted with NLO 
distributions; in this way we believe that e.g. the magnitude of K-factors more properly 
reflects the NLO subprocess terms. Also, we use throughout the NLO expression of (/i) 
with the values for the QCD scale A, flavor thresholds and number of active flavors Nif = 
N — 1 that match the definitions corresponding to heavy quark decoupling. We note 
that in ||^ the above values were taken to match the definitions for the respective parton 
distributions. However, we have explicitly verified that this amounts to a negligible change 
in the final numerical results. Note, in (2.8) we take M = m. 

At this moment there is no experimental information on the polarized photon structure 
functions AF^/^ and AFg/^, which determine the resolved 7 contributions. To estimate 



them we have used the LO maximal and minimal saturation sets of 211, as well as the sets 



of P5[, which belong to the class of the so-called asymptotic solutions. The two sets of 



21] give contributions differing little, with the maximal saturation one slightly exceeding; 



the results presented below correspond to this set. The largest resolved contributions come 



from p2 



In Figs. 31, II and HI, at = 10 and 100 GeV we present quantities related with 
the differential cross sections Ada/dpT, where pr = Pst (Fig. 1), versus xt = 2pt/\/S. 
Measurement of such cross sections at ~ 10 GeV may be carried in (a) of . Here we 

1 /2 

use the renormalization and factorization scale /i = Mf = {p\ + m^) 

In the parts (a) of Figs. 31, II and HI we present the NLO and LO (denoted by a *) 
contributions to the physical differential cross section for sets A, B and C of For set 
B we also present the contribution of subprocess (1.3) and of the resolved photon. 
In the parts (b) of the same Figs we present the asymmetries 

Ada/dpT , . 
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The unpolarized distributions are the most recent set, CTEQ5 [Q. In All the resolved 
7 contributions have been left out since they are small and what is presently known does 
not permit a completely scheme independent calculation. The errors have been estimated 
using 

- idmi 

At 75 = 10 GeV we use the conditions of [Pb = 80%, Pt = 25%, L = 2 fh-\ c-quark 
detection efficiency ec = 0.014) and unpolarized cross section a integrated over a bin of xt 
corresponding to Apr = 0.5 GeV. At = 100 GeV we use Pb = Pt = 70%, L = 100 
pb~^, ec = 0.15, for b-quark e?, = 0.05 and a integrated over a bin corresponding to Ap^ = 5 
GeV. 

Figs. 31(a) and 311(a) show that between = 10 and 100 GeV the shape of the 
LO Ad(7Lo/dpT and NLO Ada/dpr varies dramatically; this also holds for the K-factor, 
K = Ada / dpT / Ada LO / dpT ■ 

Most impotrant is the possibility to distinguish between sets A, B and C. Fig. 31(b) 
shows that at = 10 GeV near xt = 0.3 one can distinguish A and C and perhaps all A, 
B, C. Figs. 311(b) and 3111(b) show that at = 100 GeV the best range is 0.2 < xt < 0.3; 
and for Q = c one may distinguish all A, B, C, but ioi Q = b only A and C. 

In Figs. 41, II and III we present rapidity distributions. Here we use fj. = Mf = 2m. 
The presented differential cross sections are analogous to those of Figs. 31, II and III and 

, , , Ada/dY , , 



The errors have been estimated using ( |7.2| ) where now the unpolarized cross sections a are 
integrated over a bin ISY = 1. 

Fig. 41(b) shows that at = 10 GeV the region 1.25 < F < 1.5 is the best to 
distinguish set C from A or B. Fig. 411(b) shows that at = 100 GeV for c-quark, 
^ll(^) has become too small. At F ~ — 1 it seems one can distinguish all A, B, C, but 
Ada/dY is small for all sets (Fig. 411(a)). Perhaps more promising is the range < y < 1, 
where one can distinguish C from A or B. Finally Fig. 4111(b) shows that detection of 
b-quark is not useful due to large errors (e^ small). 

Figs. 51,11 and III present integrated cross sections Act and the corresponding asymme- 
tries All = Aa/a versus the cm. energy a/S*. The scale is again /j, = M = 2m. 

Comparison of Figs. 51 and 511 shows that at the two different ranges of a/S" the changes 
in the shapes and signs of Aa and All is again dramatic; clearly the same holds for the 
corresponding /^-factors, K = AaNLO / AaLo- 

In Fig. 51(b) the error (at v^=10 GeV) is estimated using again in ( [7.21) the conditions 
of [0]. Under these conditions we conclude that sets A and C can be distinguished, but not 
sets A and B or B and C. The proposed SLAG experiment p|a], which amounts to better 
conditions, and will give results at somewhat lower a/S", may distinguish also B and C. 

In Figs. 511(b) and 5111(b) the errors (at -\/S'=100 GeV) have been estimated using 
again the values of Pb, Pt, L, Ec and Eb stated after Eq. ([7.21). For c-quark, | All I are very 
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small due to relatively large unpolarized cross sections a. For the same reason, however, 
the error SA^l is not very large, so set C can be distinguished from A or B. For 6-quark, 
due to a combination of small Eb and rather small a, the error is very large and precludes 
any useful information on A^f. 

Finally, Fig. 6, for the integrated NLO cross sections Act and a and for the asymmetries 
All = ^<y/cr, shows the effect of changing the c-quark mass (part (a)) and the scales 
fi, Mf (part (b)). The results refer to set B of ||2^. E.g. regarding Act, in Fig. 6(a) we 
define 

Aa(mJ - Aa(1.5 GeV) 
Aa(1.5GeV) ' 
and in Fig. 6(b), keeping /i = Mf, we define 

Aa(/i) - A(T(2m,) 

^SC = T— 77^ ^ , (7.5) 

Acr(2mc) 

similarly for a and An. Fig. 6(a) shows that at the lower the effect of changing is 
more pronounced. 



VIII. COMPARISON WITH OTHER PUBLICATIONS 

Figs. 311(b) and 3111(b) show that at small xt, All{pt) is small; the same holds for 
of Fig. 411(b). This may lead one to conclude that HERA is rather useless in 
specifying Ag |Q. However, it may not be so. On the basis of Figs. 311(b) and Ill(b), 
reconstruct events and select only those with, say, xt > 0.2, i.e. carry integrations of 



Ada/dpT over some cut phase space. This may well enhance the resulting All 1131 •• 
course, an estimate of the corresponding errors is required to reach a definite conclusion. 

Finally, since we present analytic results for the unpolarized cross section as well, we 
will compare with similar results of ("soft" part, Eq. (2.24) of |Q]); here the relevant 
part is the last three lines of Eq. (4.5). Ref. 0] uses the phase space slicing method, which 
separates the soft and hard gluon parts via a cut parameter A. The formal relation with 
our approach is 

^ - A; (8.1) 



the necessary framework to relate these two methods is developed in (03 . Now, concerning 



terms involving t and u, we easily see that they are exactly the same, except that t and 
u are interchanged (c.f. our definition, Eq. (2.1) with that of [§], Eq. (2.13)). The only 
difference seems to arise from the coefficient of C(2), which is —2 in our case versus —3/2 
in p. Note, however, that our coefficient F{e) in (3.6) contains r(l + e), which upon 
expansion in powers of e, gives a term (£^/2)C(2); this accounts for the difference. Since 
F{e) appears both in our loop contribution (3.6) and in our Brems contribution (4.7), 
the overall result is unaffected. To verify our calculation we have evaluated numerically 
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the NLO MS scaling functions for the partonic 7(7 cross section, taking into account an 



additional "mass" factorization term given in eq. (6.31) of |]T2[, and compared it to the 
corresponding curves of Fig. 5 of 0. We found exact agreement. We have also explicitly 
verified that the sum of our non-Abelian loop contributions and the Brems ones, that are 
proportional to the Born contribution, equal analytically the corresponding "virtual+soft" 
expression presented in H]. 



IX. CONCLUSIONS 

In this paper we have presented the complete analytic results for the heavy flavor 
photoproduction for both, longitudinally polarized and unpolarized initial particles, in a 
closed form. These include the NLO contributions of the hard Brems due to the relevant 
partonic subprocesses ( |1.2|) and (|1.3|) that are presented for the first time in analytic form. 
We have computed numerically various total and differential cross sections for the energy 
ranges of CERN, SLAG and HERA. We have discussed the possibilities to differentiate 
between various scenarios for the polarized gluon distribution Agf and have once more 
emphasized the way to enhance the asymmetries for HERA energies by measuring the 
differential cross sections with the help of certain acceptance cuts (see also our earlier Ref. 



(a) of on this subject). 
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APPENDIX A 

Here we present the coefficients of the loop contributions. In the following [A]y4j, i=l,3, 
are given in App. B of 0. For [^]daa-e/ dv given in Eq. (3.2): 

Ail = AAi; Aia = -4[2(7sVt^ + 8s/t + 6)mVM + llsV^M + 24s/M + 26t/M + 

I2t^/su + st/uT - 2t^/sT]m^/T 
Ais = Ayl3/2; ^A^ = A[{2u/t-s/u)m'^/t-2s/t + 2u/s]m'^/T 
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Ai = Ai] A2^A[{2As/t-2s'^/tT + 2s/T+12t/T)m'^/tu+{s/t + 6t/s + t/T+ll) 

mh/uT - 2i7r2] 

ia = A3/2; Ai^-A[Am^s/ue -{2s''/e -s/t + t/T-2)m''/u-t/T] (Al) 
For [^dd-g/dv and [^dag/dv given in Eqs. (3.5) and (3.7): 
AA\ = m^s/e 

= A[2m^s'^/Tut^ - {Ssu/t^ + 9m/T -8-t/T- Sf/sT - st/T"^ + 2i2-u/sr2)mV'"] 
= A[m^s/Tut-2{s/t-l-t^/sT)m'^/t\ 

= -4[(16sT/t=^ + 2sV^^ + 8-M/^ + t/T)mVMT + 2(s/'u-2s/r)mVi] 
>1'3 = 4[(4sr/t2 + s/t + 2)m7'«T + 2m2sV'«t^ + 3] (A2) 

For [ls\dah/dv and [A]d(7/,/dv given in Eqs. (3.9) and (3.10): 

AE'i = 2[A{h/u + l/t)m^/t + {h/u + 7/t + Au/t^)m^ + {A/u + ?,/t){t^ + u^)/s\ 

AB'2 = -{15s^/tu-62)m^s/tu- A{s^/tu + 8)m^ /tu + 13/ A{s^/tu- 2) 

AB's = 2[(3/M + 4/t)(M-t)mVt + (l/s-lA)(t^ + M^)/M] 

ABU = 4[2(2/m + 3/t + u/t'^)m'^ - (tV^ - t + 9w + 3u^/t - 6tu/T - t^T^u + t^/T'^)/s\ 

AB'5 = 8{s'^/tu-3)m^s/tu 

AB'e = 2[2{1/T + l/U)m'^/s + 2{t^/u'^ + u^/t'^ + 2t/u + 2u/t)/s-t/uU -u/tT - 
{t^/u - 2tu + u^/t) /sTU]m^ 

B\ = 2[A{l/t-3/u)myt- {l-u/t){l/u + 2/t)m^ + At/u + 2 + 3u/t] 

B'2 = [4(29/-^^ - 6/tu + 29/t^)m^ - 2{l/u^ - l2/tu + l/t^)m^s - I3t/u + 4 - 13^^1/4 

S'a = 2[A{l/t-3/u)rn^/t+{2u/t^ -l/t-h/u)m^ -2s/u + u/t] 

B\ = -A[2{A/u+l/t)m^s/t + At/u-3s/t-2e/Tu-t/T-2t^/T'^u] (A3) 

S's = -16(l/i^ + l/'u^)m^ 

S'g = -A[{t/Uu^ + u/Tt^)rri^ - {bs/tu -l/U -l/T- t^/TUu - u^/TUt)m^ + 2tu/TU] 



APPENDIX B 

In this Appendix we list the coefficients of the Brems contributions. For AC^^ given in 
Eq. (4.3) and Ada^^/dvdw given in Eq. (4.5), the coefficients Ae^ and A/j are 

Aei = 16[A{s/s2U — 1/u + l/t)'m'^ /t + 3s/s2U — 2/u — s/s2t — S2/tu — u/s2t\ 
Ae2 = S[S{-s/u-l + s/t)m'^/s2t + 2{2s/u-2 + As/t-3ss2/t{s + u))m^/s2- 
As^ / S2U + As/u — 2s/s2 — u/s2 + u/t + 2ss2/tu + 3u/ {s + u)] 
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Aes = 0; /^.e^ ^ %[A{s/u-\- s/ S2 - t/ S2)rn^ /t -\- ss2/tu -\- su/ S2t -2u/ S2 -2] 

Ae5 = d>[9>{t/u+l)m^/s2t-2{s^/s2U + s/u + At/s2-2s2/u)m^/t + st/s2U + 

4s/s2 - 2 + su/s2t + 2u/t\ 
Aee = A[8{s/u — S2/t + u/t)rn^/s2 + 2{2st/s2U + s/s2 — S2/U + 2 + 2ss2/tu — 

2S^/S2t)m'^ + + 1*^)752^] 

Ae7 = 4[8(s2/m- l)mV^ + 2(25/1^ + 3 -S2/s-2s2Ati + 2s/i)m^s/(s + ii) - 

ss2(s^ + s\)/t{s + 1^)1*]; Aeg = 16(s2/ii - u/s2)/t 

Aeg = 4[8(t/M + s/t)m^/s2-2(2s^/s2M-2s/M + 2s2/M-2 + (st/s2 + 2su/s2- 

S2)/(s + M))m2 - st{2s/s2U - 2/u + (52/^ + u/s2)/{s + li))] 
Aeio = I6{t/u — u/t) / S2; Aeii = —16m'^ s / u 

A/i = 16[4(s/K-l)mVs2t + (3/s2 + lA)sH; A/2 = 32(s + M)/is2^i 

A/3 = A/2; AU = -32[i2/s2 + 2/t~s/s2u)myt+l/t+l/s2] 

A/5 = 8[2(s/s2M + s/tM - 2u/t'^ + s/t(s + m) - 4ss2A^(s + u))m'^ - / S2U + 

2s /u - 2s/ 32 + 2- 2u/s2 + ss2/tu + S2A + (^^ + 2s2u/t)/{s + -u)] 
A/e = 4[2{2s^/s2U + s/u + S2/u + 3s/s2 + {s + 2su/s2-3s2-4:s{sl + 

u^) jtu) /{s + u))m^ + s^/s2 + 2su/s2 + 2u^ / S2 - 2s^/t - Asu/t + 

2S2A - 2^7^ + {si + -u^/ls + '")] 
A/7 = 0; A/8 = 

A/g = 4[2(2s/m - 2S2AM + S2A - (SS2/M + 2sM/t - S2M/t)/(s + M))m^ + S2 - 

2t - sl/t - v?lt + (s2 + + m)] 

A/10 = -8(2s + t)m'; A/n = (Bl) 

For G'"'^ given in Eq. (4.3) and da^^jdvdw given in Eq. (4.5), the coefficients and /j are 

ei = 16[4(l/s2'u + l/s2t-lA^)m^-3s/s2'u + 2/'u-l/s2-S2Aw-lA] 

62 = 8[-16(l/t+l/M)mVs2t-8(3t/s2M + 2/M + 2/s2- l/(s + M))m7t + 

2(sVs2M - 8s/m + 4s/s2 - 1 + 2m/s2 - 3S2A + u/t)m^/{s + u) + u/s2 - 
Ast/s2U - 2s /s2 + (2s2 -2u + 2ssl/tu + S2u/t)/{s + u)] 

63 = — 16m^ 

64 = -8[8m^(l/ii+ I/S2) - 2m7s/w- s/s2 - 2s2/ii- 2ii/s2) + (s2 + 'u)(2 + 

2ii/s2 + 2s/s2 - s/u + s/s2)]/t 

65 = 8[8{s/s2 — l)rn^/tu — 2m,'^{s/s2U + 2/u — 5/s2 + 3/t — u/ts2) + /s2U — s/u — 

S/S2 + 4 + S'u/tS2 - 2S2/t + 2s'^/tS2] 

66 = A[lQ>m^{l/su - I/S2U - l/ts2) + 8{su - {s - S2Y - {s2 - uf)m^/ts2U + 2m'^(t/u + 

S/S2 - 1 + s/t) - S{S^ + U^)/S2t\ 

67 = 4[16(s2/ii - l)m^/st + 8(s^ + SS2 + S2)?^V^^(« + ^1) + 2m^(3s/t + {ss2/u - S2 - 
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2ss2/t)/{s + u)) — ss2{s^ + sl)/tu{s + ul)]; es = 16(s2 + uf /ts2U 

eg = 4[16(s/s2 — 2 + S2/ s — u/s)mP/tu + 8{s2/u — 3st/s2U — 2 + u/ S2)m^/ {s -\-u) -\- 2m^ x 

((s - S2){l/u + I/S2) - 1 + 2(s2 - s)/{s + u)) + {S2/U - 2 + ii/s2 - 2st/s2u)st/{s + u)] 
eio = -16{l/t-l/u){s/s2-l); en^^^m^/u 
fi = -16[Sm'^/ts2U + 2{2/u + s/tu)m'^/s2 + s/u{l/t + 3/s2) - 2ss2/t{s + u)u] 

f2 = 32(l/i+l/s2)/li; /3 = /2 

/4 = -32[4m^/is2'u + 2(l/s2 + l/t)mV'"- l/s2 + S2A(s + 'u)] 

/5 = -8[8m\l/s2U + 1/tu + {1/t - 2s2/t^)/{s + u)) - Am^l/t + 2s2/t^ - {st/s2 + 

ss2/t)/{s + u)u) -t/u + 2t/s2 - sl/tu + {sH/s2U + 2u - 2s2u/t) /{s + u)] 
/e = A[Sm^{s^ /tS2U - 1/S2 - A/t - {s/t - S2s/ut - S2u)/{s + u)) + A{s^/s2 + su/s2 - 

2S2 + S-U- S2{S + S2)/{S + U))m^/t - (1 - {si + U^)/{S + Ul)t){s + S2 + uf/s2] 

/7 = -32(s + M)W/t2; /s = 32(s + M)mV^ 

/g = A[Sm^{s^/tu + S2/u + S2/t)+Am^{ss2/t + S2-2s) + {s2/t-2){sl + u^ -st- 

ut)]/{s + -u); /lo = lQm^{2m^ + -u); /n - -Sm^^ (B2) 

Now we shall write down the coefficients [AJcj and [A]/j for the subprocess 75 — > QQ?- 
For Qi we have: 

Aei = -8(2m^-^)/^^ Aeg = 0, Ae4 = -8*s/i, Aeg = 
A/4 = 8(2m2 + ^)/^^ A/5 = 4(2(s2 + s)m2-iii)/i2 

A/6 = -2(2m' + 2s2-i), A/y = 0, A/g = 0, A/10 = 2(2s + ^)m^ A/n = 
ei = 8(2m^ + i)/i^ 63 = 4m^ 64 = 8(2w + eg = -16/i 

/4 = /s = -4(4m^ + 4m2s2 + ii^)A^ 

/e = 2(4(s2 + u)m^ + 2(2s + ?>t)s2 + 8m^ - 4.sl - 2s^ - Ast - 3t^) /t 
fr = 8{s2-tfmyf, fs^-8u{s2-t)myt 

/lo = -4(m' + ^^)m^ /u = 2iiW (B3) 



For 



2- 



Aei = -16/s, A/i2 = -8/s, A/13 = -SmVs, Af^ = -8{2m^ + u)/s 

A/15 = 8t/s, A/16 = 2(2m(s + t)- 4m^s + + 2ts2) /s 

A/17 = 2m2s, A/2o = 

ei = 16/s, /i2 = 8/s, /i3 = 8mVs, /14 = 8(2771^ + 

/i5 = 8{2m^ + m) /s, /i6 = -2(4m^s - 4s^ + 6S2S + 4^2^ - 3s^ - Ast - 2t^) /s 

/i7 = 2m's, /2o = -16/5 (B4) 

And, finally, for Q^: 

Aei = 8{s2-2s + t)/st, Ae4 = -8, Aeg = 
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A/4 = -8(s2 - t) /st, A/s = 4(2(s2 - t)m'^ + si - S2t - ut) / st 

A/e = -2{2{2s + t)m'^ -2s2t + t^)u/{s + t)s, A f 12 ^ S{s2 - s) / st 

A/14 = 4(2(s2 - 3s)m^ - 2su + 2s2U + ^52) /st, A/15 = 8 

A/16 = -2{2{s + 2t)s2-2m'^s- s^ -4.st-2t^)u/{s + t)t 

A/18 = -2{2{{2s + t)s2 + 2s^ + 2st)m'^ - {2u + t)s2t)/{s + t)s 

A/19 = -2{2{s + t)s2 - 2n?s - + 2ut)s2/{s + t)t, A/20 = 0, A/21 = 

ei = -8(s2 -2s-t) /st, 64 = 8(2m^ + S2 + U + tu/(s + t)) /s, eg = -16/s 

h = -8{s2-t)/st, h = A{2{s2-t)m^ + sl-S2t-ut)/st 

fe = -2{2{As/t + t/{s + t))m'^ + 2s2 + 2u + {2s2t-4:sl-t'^)/{s + t))u/s 

fi2 = -8(s2 - s) /st, /i4 = -4(2(s2 - 3s)m2 + 2s2^i + - 2su) /st 

/i5 = 8(2m2 + M + ss2/(s + t))/t 

/16 = 2(2m2s-4s2 + 6s2S + 4s2t-3s^-4st-2t^)M/(=s + ^)^ 

/18 = -2(2(4sM/t + 2s - S2t/(s + t))m'^ + 52(^2 + « - ((s - 2s2)^ + 2s2s)/(s + t)))/s 

fi9 = 2(2(s + 2t)s2 -2mls-4s2 -s^ -2st-2t2)s2/(s + t)i 

/20 = -16/t, /2i = -16m2s/t (B5) 



APPENDIX C 



We give here the Brems integrals, F^, i = 1,...,11, appearing in Eq. (4.5) and also in 
Eqs. (5.4) - (5.6). Define, 

Pi^us2-s{t + 2m'^), P2^ s{{t + m^){us2-m^s) -m^{s + uf) (CI) 

We may now write down the integrals: 

„ 1 r„ SS<2 , o, , So , 

Fi = —, r^[Pi + — (i + 2m')ln 



2 J 



(s + uY S2 m 

[s + uy Ab2 2b2 S2 

.sHt + 2m^ m2 5 
2 So 



S2S2U 



2 



2(s + m)3 



S2 r{t + uf 2 t + u + 2m^ 3 Pf 



F| = TTTTT^. — -[- — r^-m^s ^ Pi-- 



2^2(5 + m) 4 s + M 4(s + m)2' 

^ _S2V_ c ^ Pi 2S2 

^ {s + uy' ^ {s + uy' ^ S2{s + uy ' 
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S2ut^ S2ut rn? SS2 + ut 



165*9 So S -\- t , n, 

^ {s + u){ut-2m^s) 



slut{s + u) 2S2 uH^ 



165*0 r So S -\- t , O \ /, So , ut . So 2Po -, 



s^u^ S2U rn? {s + uY 

^ 8^2(g + M) ^ 8g2(g + M) M^ '^2 ^ ^ g2 Mt-2m^g 

S2UH^ ' S2UH^ m? SS2 + "Ut SS2 + 



sluH"^ SluH"^ S2 m? SS2 + ut Si uH"^ 

Note that parts of , Fg and Ff^ proportional to 1 — e cancel out exactly in Eq. (4.5). 
The singular and finite parts of these integrals can be found in Appendix C of [Q; below 
we give the derivation of 0{e) terms. 

We use the momentum parametrizations of Appendix A of and, as in |jl2|, we denote 

cos + C sin 6*1 cos 6^2) ' (C3) 

where ^ ^ 

/ d9ism^-^'9i I d92sm-'^'92. (C4) 
J Jo Jo 

All the above integrals are proportional to 1/a^ = 1/ujIuj2 ~ 1/^2 ~ 1/(1 — w)^; since the 
2^3 particle phase space is proportional to (1 — it?)^^^^, in view of Eq. (4.6), the terms 
of 0{e) give finite contributions proportional to 5(1 — w). 

Integral F7 = /c?f2„/(p2 ■ k^'- 

This is of the type I^''^^ of 0]. The result is 

F7 = -^T^ = -^{l-e + 0{e')) (C5) 



1 + e a 



Integral Fg = JdQn/ {P2 ■ kfi.Pi ■ Pa)- 
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This is of the type 1^'^^ of [0, and determination of the 0{e) term proceeds as follows. 
First, defining 

H = A + Bcosei + Csmeicose2, (C6) 

one can show the identity 

1 1 ( 1 B Csin^icos^sl 



(l-cos^i)iJ A + B [1-cosei H (l-cos^i)// 
and by repeated application of it: 

sin 6i 
(1 -cos^i)2i7 ~ 

1 sin^^i B |. sin^i ^ i?sin6'i C(l + cos^i) cos^2-| 

171 „„„ a ~'~ 



A + 5'(l-cos^i)2 A + -cos^i) H H 



Ccos^2r(l + cos^i) 5(1 + COS .sin^' COS^ ^2 



+ 



A + 5 ^(1 -cos^i) H ' (A + 5)2^(l-cos^i; 



2 



Ssin6'i(l + cos6'i) cos^6'2 C(l + cos6'i)2 cos^ 6'2,, . , 

+ jj ]} (C8) 

Note that the fifth term vanishes due to the integration over 62. Also, terms with only 
H in the denominators are finite, consequently have no poles and cannot produce finite 
contributions from their 0{e) terms. Thus, we are left with the terms 1,2 and 7. After 
integrating them in n-dimensions, summing and keeping the relevant order e terms, we 
arrive to the following result: 

Integral Fu = Jdnn/{P2 ■ kfiPi -PiY- 

Proceeding as before, the term of 0{6) is provided by 

Finally, we give the Brems integrals, Fi, i = 12, 21, appearing in Eqs. (5.4) and (5.6). 
With 

Ps = 2m^s + u{s- S2), Z = Am^s{s + t) - sit (Cll) 



we have: 



F12 = -^[Piit^u) + ^P^Fu] 

y ^ 
1 2 s* 

Fxz = --[^P3 + Piit^u)F,^] 

yZ gjy^Z 
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2S2 ( S2{s - S2) + 2m^s + S2y\ 
^S2{s - S2) + 2w?s — S2y J 




Fi4 = 



S2{s + t)' 

^16 = , ^16 = In— + ln 

852(^ + t) 852(. + t) ^2 «^ S2Ps 

^17 = „ „2„,2 ' ^17 = „ o .2 i^"Z:2 



2 



S2S^u^ ' S2s'^u^ m? {s + 1)^ m?s{s + 

S2\/—tz — 1vn?s{s + — S2\/—tz J 

s _ ^^2 c _ 8'S'2 , S'S'2 

-^19 - ~T27' -^19-727^^ 



S2S S2S SS2 + Mi 

p. _ ^2^ _ Pl(t^M) 
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FIGURE CAPTIONS 

Fig. 1. LO and loop graphs. In the loop graphs pi ^ p2 crossed ones are not shown. Note 
that graph (i), representing gluon, quark and ghost loop, does not contribute here. 

Fig. 2. A) Gluon Brems graphs; pi <->■ p2 crossed ones are not shown. B) Graphs of the 
subprocess •yq QQq. 

Fig. 3. Quantities related with the pr distributions versus xt = 2pT/\fS: Parts (a): Po- 
larized differential cross sections; the LO (Born) ones are indicated by *. Parts (b): 
Asymmetries for sets A, B and C. 31: Q = c, = 10 GeV. 311: Q = c, = 100 
GeV. 3III: Q = 6, 75 = 100 GeV. 

Fig. 4. Quantities related with the rapidity Y distributions: Parts (a) and (b), as well as 41, 
411 and 4III as in Fig. 3. 

Fig. 5. Quantities related with the integrated cross sections for 7p ^ Q + X: (a) Factors 
K = Aa/Aa^o (b) Asymmetries. 

Fig. 6. At cm. energies = 10 and 100 GeV, for integrated cross sections and with solid 
lines for Aa, dashed for a and dotted for the asymmetry All. a) The ratio Rm (see 
end of Sect. VI) with m = mc', b) The fractional variation Rgc with the scale /i = Mf 
and with respect to /i = Mf = 3 GeV. For both a) and b) the lines specified by + 
refer to the corresponding quantities for -s/S* = 100 GeV. 
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